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Abstract 


This  report  describes  the  application  of  the  chimera  numerical  technique  to 
a  multi-body  segmented  projectile  configuration  system  of  interest  to  the 
U.S.  Army.  Computations  were  performed  at  a  supersonic  speed  on  this 
configuration  which  consists  of  an  ogive-cylinder  projectile  with  a  peg¬ 
shaped  trailing  segment.  The  computed  results  show  the  qualitative 
features  of  the  wake  flow  field  for  the  projectile  with  the  segment  in  three 
different  positions:  centered,  offset,  and  angled.  The  segment  in  the  offset 
position  has  a  strong  effect  on  the  flow  field  in  the  aft  region  of  the 
projectile,  thus  affecting  the  aerodyamic  coeffiecients  of  the  projectile.  The 
force  and  moment  coefficients  of  Ae  segment  are  also  significantly  affected 
by  the  orientation  of  the  segment. 
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COMPUTATIONAL  MODELING  OF  A  SEGMENTED  PROJECTILE 


1.  INTRODUCTION 

Aerodynamic  forces  and  moments  are  critical  design  parameters  used  in  the  design  of  artillery 
shells  and  bodies  flying  in  relative  motion  to  each  other.  The  advancement  of  computational  fluid 
dynamics  (CFD)  is  beginning  to  have  a  major  impact  on  projectile  design  and  development. [1-4] 
Improved  computer  technology  and  state-of-the-art  numerical  procedures  enable  scientists  to 
develop  solutions  to  complex,  three-dimensional  (3-D)  problems  associated  with  projectile  and 
missile  aerodynamics.  The  research  effort  has  focused  on  the  development  and  application  of  a 
versatile  overset  grid  numerical  technique  to  solve  geometrically  complex  single-body  as  well  as 
multi-body  aerodynamic  problems.  This  numerical  capability  has  been  used  successfully  to 
determine  the  aerodynamics  on  a  number  of  multi-body  projectile  configurations  [4-6]  at  transonic 
and  supersonic  speeds.  Earlier  applications  involved  axisymmetric  flow  computations.  [4,5] 
Recently,  this  technique  has  been  used  to  investigate  the  submunition  dispersal  from  an  Army 
tactical  missile  system  [6]  involving  3-D  flow  computations.  This  report  describes  the  application 
of  the  advanced  numerical  technique  to  a  multi-body  segmented  projectile  configuration  of  interest 
to  the  U.S.  Army.  Figure  1  shows  a  computational  model  for  this  system,  with  several  trailing 
segments  centered  along  the  line  of  symmetry  in  the  wake  of  the  projectile.  This  multi-body 
problem  involves  3-D  flow  computations  of  trailing  segments  flying  in  the  wake  of  a  parent 
projectile.  However,  the  scope  of  this  study  is  limited  to  one  trailing  segment.  The  particular 
problem  here  is  to  determine  the  resulting  aerodynamic  interference  with  the  segment  in  centered, 
offset,  and  angled  positions. 


Figure  1.  Computational  Model  for  a  Segmented  Projectile. 


The  complexity  and  uniqueness  of  this  type  of  multi-body  problem  result  from  the 
aerodynamic  interference  of  the  individual  components,  which  include  3-D  shock-shock 
interactions,  shock-boundary  layer  interactions,  and  highly  viscous-dominated  separated  wake 
flow  regions.  The  overset  grid  technique  [7-9],  which  is  ideally  suited  to  this  problem,  involves 
generating  numerical  grids  about  each  body  component  and  then  oversetting  them  onto  a  base  grid 
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to  form  the  complete  model.  With  this  composite  overset  grid  approach,  it  is  possible  to 
determine  the  3-D  interacting  flow  field  of  the  multi-body  system  and  the  associated 
aerodynamic  forces  and  moments  at  different  positions  and  orientations  without  the  need  for 
costly  regridding.  The  solution  procedure  of  the  developed  technique  is  to  compute  the 
interference  flow  field  at  multiple  locations  until  final  converged  solutions  are  obtained  and  then 
to  integrate  the  pressure  and  viscous  forces  to  obtain  the  total  forces  and  moments. 

A  description  of  the  computational  algorithm  and  the  chimera  technique  follows.  The  next 
section  describes  the  model  geometry  and  various  computational  grids  used  in  the  numerical 
computations.  Steady  state  computational  results  are  presented  for  various  orientations  and 
locations  of  the  segment  in  the  projectile’s  wake. 


2.  SOLUTION  TECHNIQUE 

The  complete  set  of  time-dependent,  Reynolds-averaged,  thin  layer  Navier-Stokes 
equations  is  solved  numerically  to  obtain  a  solution  to  this  problem.  The  numerical  technique 
used  is  an  implicit,  finite  difference  scheme.  Time-accurate  calculations  are  made  to  numerically 
simulate  the  submunition  in  three  possible  positions  relative  to  the  main  projectile. 


2.1  Governing  Equations 


The  complete  set  of  3-D,  time-dependent,  generalized  geometry,  Reynolds-averaged,  thin 
layer  Navier-Stokes  equations  for  general  spatial  coordinates  T],  and  ^  can  be  written  as  follows 


d^q  +  d^F  +  dj^G  +  d^H  =  Re'^d^S , 

where 

^  =  4(x,  y,  z,  t)  -  longitudinal  coordinate; 

Tj  =  T|(x,  y,  z,  t)  -  circumferential  coordinate; 

^  =  ^(x,  y,  z,  t)  -  nearly  normal  coordinate; 

T  =  t  -  time 


(1) 


In  Equation  1,  q  contains  the  dependent  variables:  density,  three  velocity  components,  and  energy. 
The  thin  layer  approximation  is  used  here,  and  the  viscous  terms  involving  velocity  gradients  in 
both  the  longitudinal  and  circumferential  directions  are  neglected.  The  viscous  terms  are  retained  in 
the  normal  direction,  ^ ,  and  are  collected  into  the  vector  S.  These  viscous  terms  are  used  every¬ 
where.  In  the  wake  or  the  base  region,  similar  viscous  terms  are  also  added  in  the  streamwise 
direction,  An  implicit,  approximately  factored  scheme  is  used  to  solve  these  equations. 
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2.2  Numerical  Algorithm 


The  implicit,  approximately  factored  scheme  for  the  thin  layer  Navier-Stokes  equations 
using  central  differencing  in  the  Tj  and  ^  directions  and  “upwinding”  in  ^  is  written  in  the 

following  form  [11]: 


1  + 


ii,hdl{A-^ )”  +  ithd^C'^  -  ibh  Re-^5i;J-^MV  -  ^ 


I  +  ii)h8^  (^  >  1 7] 


=  (G”  -G„) 

Re-% (5”  -L)}- ibDe[Q'’ " \ 


(2) 


where  h  =  At  and  the  free-stream  base  solution  is  used.  Here,  5  is  typically  a  three-point  second 
order  accurate  central  difference  operator,  5  is  a  midpoint  operator  used  with  the  viscous  terms, 
and  the  operators  and  are  backward  and  forward  three-point  difference  operators.  The 
flux  F  has  been  eigensplit,  and  the  matrices  A,  B,  C,  and  M  result  from  local  linearization  of  the 
fluxes  about  the  previous  time  level.  Here,  J  denotes  the  Jacobian  of  the  coordinate 
transformation.  Dissipation  operators  De  and  Dj  are  used  in  the  central  space  differencing 
directions.  The  smoothing  terms  used  in  the  present  study  are  of  the  form 


and 

where 


-7 


\s2Sp{B)pd  +  e,5  I,/, 


£i,|,  =  (A/y-'[£j«p(5);85  +  2.5£4«p(B)5]  I,  J, 


and  p(B)  is  the  true  spectral  radius  of  B.  The  idea  here  is  that  the  fourth  difference  will  be  tuned 
down  near  shocks  (e.g.,  as  p  gets  large,  the  weight  on  the  fourth  difference  drops  down  while  the 
second  difference  tunes  up). 
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2.3  Chimera  Composite  Grid  Scheme 

The  chimera  overset  grid  scheme  [7-9]  is  a  domain  decomposition  approach  in  which  a 
configuration  is  meshed  by  using  a  collection  of  overset  grids.  It  allows  each  component  of  the 
configuration  to  be  gridded  separately  and  overset  into  a  main  grid.  Overset  grids  are  not  required 
to  join  in  any  special  way.  Usually,  a  major  grid  covers  the  entire  domain  or  a  grid  is  generated 
about  a  dominant  body.  Minor  grids  are  generated  about  the  rest  of  the  other  bodies.  Because 
each  component  grid  is  generated  independently,  portions  of  one  grid  may  lie  within  the  solid 
boundary  contained  within  another  grid.  Such  points  lie  outside  the  computational  domain  and 
are  excluded  from  the  solution  process. 

In  the  segmented  projectile  study,  the  grids  around  the  projectile  body  comprise  the  major 
grid,  while  the  grid  around  the  segment  is  the  minor  grid.  The  minor  grid  is  completely 
overlapped  by  the  major  grid^  thus,  its  outer  boundary  can  obtain  information  by  interpolation 
from  the  m^or  grid.  Similar  data  transfer  or  commumcation  is  needed  from  the  minor  grid  to  the 
major  grid.  However,  a  natural  outer  boundary  that  overlaps  the  two  grids  does  not  exist.  The 
chimera  technique  creates  an  artificial  boundary  (also  known  as  a  hole  boundary )  between  grids, 
which  provides  the  required  path  for  information  transfer  from  the  segment  to  the  projectile  grid. 
The  resulting  hole  region  is  excluded  from  the  flow-field  solution  in  the  projectile  grid.  Equation 
2  has  been  modified  for  chimera  overset  grids  by  the  introduction  of  the  flag  ib  to  achieve  just 
that.  This  ib  array  accommodates  the  possibility  of  having  arbitrary  holes  in  the  grid.  The  ib 
array  is  defined  so  that  ib  =  1  at  normal  grid  points  and  ib  =  0  at  hole  points.  Thus,  when  ib  =  1, 
Equation  2  becomes  the  standard  scheme.  When  ib  =  0,  however,  the  algorithm  reduces  to 

=  0  ox  =  0”,  leaving  Q  unchanged  at  hole  points.  The  set  of  grid  points  that  forms 
the  border  between  the  hole  points  and  the  normal  field  points  are  called  inter-grid  boundary 
points.  These  points  are  updated  by  interpolating  the  solution  from  the  overset  grid  that  created 
the  hole.  Values  of  the  ib  array  and  the  interpolation  coefficients  needed  for  this  update  are 
provided  by  a  separate  algorithm.  [7] 

A  major  part  of  the  chimera  overset  grid  approach  is  the  information  transfer  from  one  grid 
into  another  by  means  of  the  inter-grid  boimdary  points.  Again,  this  set  of  points  defines  the 
hole  boundaries  and  outer  boundaries  of  the  minor  grids.  These  points  depend  on  the  solutions  in 
the  overlapping  regions.  In  the  present  work,  the  PEGSUS  code  [12]  has  been  used  to  establish 
the  linkages  between  the  various  grids  that  are  required  by  the  flow  solver  or  aerodynamics  code 
described  earlier.  These  include  the  determination  of  the  interpolation  coefficients  and  the 
establishing  chimera  logic  for  bodies  making  holes  in  overlapping  grids. 
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2.4  Boundary  Conditions 


For  simplicity,  most  of  the  boundary  conditions  have  been  imposed  explicitly.[l]  An 
isothermal  wall  boundary  condition  is  used  on  the  body  sxuface,  and  the  no-slip  boundary 
condition  is  used  at  the  wall.  The  pressure  at  the  wall  is  calculated  by  solving  a  combined 
momentum  equation.  Free-stream  boundary  conditions  are  used  at  the  inflow  boundary  as  well 
as  at  the  outer  boundary.  A  symmetry  boundary  condition  is  imposed  at  the  circumferential 
edges  of  the  grid,  while  a  simple  extrapolation  is  used  at  the  down-stream  boundary.  A 
combination  of  symmetry  and  extrapolation  boundary  conditions  is  used  at  the  center  line  (axis). 
For  cases  where  the  free-stream  flow  is  supersonic,  a  nonreflection  boundary  condition  is  used  at 
the  outer  boundary.  A  no-slip  wall  boundary  condition  is  used  on  the  body  surface  of  the  trailing 
segment.  Boundary  condition  is  not  applied  at  the  outer  boundary  of  the  trailing  segment; 
instead,  it  is  updated  through  the  chimera  interpolation  procedure. 


3.  MODEL  GEOMETRY  AND  COMPUTATIONAL  GRID 

The  projectile  model  consists  of  an  ogive-cylinder,  6  calibers  long,  with  a  peg-shaped 
trailing  segment.  The  segment  is  1.5  calibers  long  and  is  located  1  caliber  behind  the  projectile. 
Figure  2  shows  a  cross-sectional  view  of  the  computational  grid  for  the  multi-body  configuration. 
The  entire  grid  system  consists  of  approximately  2  million  grid  points  and  is  split  into  five 
sections:  a  small  grid  in  front  of  the  projectile  body,  two  grid  zones  along  the  projectile  body,  a 
large  grid  in  the  wake  region  of  the  projectile,  and  finally,  a  grid  around  the  trailing  segment.  Each 
grid  section  was  obtained  separately  and  then  combined/overset  to  provide  the  full  grid.  The 
total  number  of  points  in  the  longitudinal,  circumferential,  and  radial  directions  for  the  projectile 
grids  is  156x75x70.  The  projectile  wake  or  the  base  region  grid  is  94x75x129  and  the  segment  grid 
is  1 1 1x75x25.  The  projectile  grids  extend  2  calibers  in  tlie  nomal  direction  from  the  line  of 
symmetry,  while  the  segment  grid  extends  only  0.1  caliber  from  the  segment  surface  and  is 
entirely  embedded  in  the  base  region  grid.  A  body-conforming  grid  was  obtained  for  the  segment 
and  then  overset  to  form  the  composite  mesh  shown  in  Figure  2.  This  case  corresponds  to  the 
centered  position,  where  the  segment  is  centered  on  the  line  of  symmetiy,  behind  the  projectile. 
The  major  grid  comesponds  to  the  grids  for  the  projectile,  including  the  wake,  which  are  easily 
generated  independently  of  the  minor  grid  (grid  for  the  segment).  In  this  and  all  subsequent 
figures,  a  portion  of  the  segment  surface  has  been  cut  away,  allowing  a  view  of  grid  and/or  flow 
field  in  the  cavity. 
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Figure  2.  Computational  Grid  System. 


As  stated  earlier,  the  chimera  technique  allows  individual  grids  to  be  generated  with  any  grid 
topology,  thus  making  the  grid  generation  process  easier.  For  this  study,  the  segment  was  initially 
positioned  in  the  wake  of  the  projectile,  centered  along  the  longitudinal  line  of  symmetry.  A 
geometric  transformation  was  used  to  place  it  in  both  an  offset  and  an  angled  position  for  additional 
configurations.  Figure  3  shows  an  expanded  view  of  the  grid  for  each  segment,  embedded  in  the 
projectile  wake.  Again,  there  was  no  need  to  re-grid  to  generate  the  minor  grid;  a  geometric 
transformation  was  used  to  place  the  original  (centered)  grid  in  the  desired  position.  As  part  of  the 
chimera  blanking  procedure,  the  segment  cuts  a  hole  in  the  base  region  grid  of  the  projectile  (see 
Figure  4).  All  grid  points  in  the  base  region  grid  that  lie  within  this  hole  are  blanked  and  are  excluded 
from  the  flow  field  solution  process.  Also,  all  points  on  this  hole  boundary  in  the  base  region  grid 
are  updated  through  interpolation  of  the  solution  in  the  segment  grid.  The  grid  blanking  (hole 
cutting)  shown  in  Figure  4  corresponds  to  the  segment  in  the  centered  position.  Similar  hole  cutting 
is  done  for  other  positions  of  the  segment  in  the  wake  of  the  parent  projectile. 


4.  RESULTS 

In  this  study,  3-D  steady  state  numerical  computations  were  performed  for  the  projectile 
alone  and  with  a  trailing  segment  in  the  three  positions  of  interest  described  earlier.  All 
computations  were  run  at  Moo  =  2.5  and  a  =  0°,  and  atmospheric  flight  conditions  were  used. 
The  calculations  involving  the  trailing  segment  required  approximately  45  million  words  of 
memory,  and  each  case  used  an  average  of  30  hours  of  computer  time  on  the  Cray  C-90 
supercomputer  at  the  Aeronautical  Systems  Center  major  shared  resource  center  (MSRC).  The 
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calculation  for  the  projectile  by  itself  used  fewer  resources  and  was  done  at  the  U.S.  Army 
Research  Laboratoiy  MSRC  on  a  Silicon  Graphics  platform. 
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Figure  4.  Background  Grid  Showing  Chimera  Blankint 
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Results  are  shown  first  for  the  case  when  there  is  no  trailing  segment.  Figures  5  and  6  show 
Mach  contours  and  pressure  contours,  respectively.  Here,  blue  represents  low  values  and  red 
represents  high  values  of  Mach  and  pressure  contours.  These  figures  clearly  show  the  expected 
flow  features.  Figure  5  shows  an  oblique  shock  wave  emanating  from  the  nose  of  the  projectile. 
The  flow  expands  at  the  base  corner,  which  is  followed  by  a  recompression  shock  down  stream 
from  the  base.  Figure  6  shows  high  pressures  in  the  nose  region  and  low  pressures  near  the  base. 
As  expected,  the  flow  field  is  quite  symmetrical. 


Figure  5.  Mach  Contours  for  a  Projectile  Without  Segments. 


Figured.  Pressure  Contours  for  a  Projectile  Without  Segments. 


Computed  results  obtained  for  the  projectile  and  the  segment  are  presented  next.  The 
segment  is  placed  in  the  wake  of  the  projectile  in  the  centered,  offset  (from  center  line  of 


symmetry),  and  angled  (segment  at  angle  of  attack)  positions.  Figure  7  shows  Mach  contours  for 
each  of  the  three  segment  orientations.  Similarly,  Figure  8  shows  pressure  contours  for  each  of 
the  three  segment  orientations.  As  seen  in  these  figures,  the  flow  field  in  the  wake  region  has 
changed  dramatically  because  of  the  presence  of  the  trailing  segment.  As  expected,  the  flow  field 
in  the  wake  is  still  symmetrical  when  the  segment  is  in  the  centered  position.  The  asymmetry  in 


Figure  7.  Mach  Contours  for  a  Segment  in  Three  Positions:  Centered.  Offset,  and  Angled  (top 
to  bottom). 
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Figures.  Pressure  Contours  for  a  Segment  in  Three  Positions:  Centered.  Offset,  and  Angled 
(top  to  bottom). 


the  wake  region  flow  can  be  clearly  seen  when  the  segment  is  in  the  offset  and  angled  positions. 
There  has  been  a  significant  increase  in  the  width  of  wake  region  flow  for  the  case  of  the  segment 
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in  the  offset  position  (see  Figure  7).  In  general,  the  pressure  is  higher  in  front  of  the  segment  and 
lower  in  the  base  region  of  the  segment  itself.  Of  the  three  orientations  of  the  segment,  the  offset 
case  seems  to  have  the  largest  effect  in  the  flow  region  between  the  projectile  and  the  segment. 

As  seen  in  Figure  7,  the  flow  field  in  this  region  is  quite  asymmetrical.  In  addition,  the  segment  in 
the  offset  position  seems  to  have  a  strong  effect  on  the  flow  field  in  the  aft  region  of  the 
projectile  upstream  from  the  base  comer  (especially  in  the  bottom  side). 

Computed  surface  pressures  on  the  trailing  segment  are  quite  interesting.  Figure  9  shows  a 
front  view  of  the  segments,  including  the  peg  and  cylinder  of  the  segment  in  each  of  the  three 
positions.  When  the  segment  is  centered  on  the  line  of  symmetry,  there  is  an  area  of  high 
pressure  along  the  inner  and  outer  edges  of  the  face  of  the  segment  behind  the  peg.  When  the 
segment  is  offset,  the  pressure  along  the  face  of  the  segment  cylinder  is  shown  to  increase 
dramatically,  especially  along  the  upper  surface  of  the  segment  peg.  Finally,  when  the  segment  is 
in  the  angled  position,  the  pressure  is  only  slightly  increased  on  the  face  of  the  cylinder,  with  a 
much  more  dramatic  change  along  the  cylinder  of  the  segment  itself  Figure  10  shows  a  rear  view 
of  the  surface  pressure  in  the  cavity  of  the  segment  in  each  of  the  three  orientations  of  the 
segment.  In  general,  the  pressure  in  the  cavity  is  low  for  all  three  segment  orientations,  being  just 
slightly  higher  when  the  segment  is  in  the  offset  position. 

Aerodynamic  forces  and  moments  have  been  obtained  from  the  computed  solutions.  Tables 
1  and  2  show  some  of  the  computed  aerodynamic  force  and  moment  data  for  the  segmented 
projectile  system,  based  on  the  position  of  the  segment.  Both  tables  show  the  axial  force  (Ca)» 
normal  force  (Cn),  and  pitching  moment  coefficients(CiTip).  The  presence  of  the  segment  in  the 
wake  of  the  projectile  is  expected  to  change  these  coefficients  both  for  the  projectile  and  the 
segment.  Table  1  shows  these  coefficients  for  the  projectile.  The  axial  force  is  only  slightly 
affected  by  the  segment  orientation;  all  three  values  are  approximately  10%  lower  than  the 
computed  value  (0.32)  for  the  projectile-alone  case.  For  the  segment  in  the  centered  position,  the 
flow  field  as  explained  earlier  is  symmetric,  and  therefore,  the  normal  force  and  pitching  moment 
coefficients  are  zero.  The  effect  of  the  angled  segment  on  the  projectile  force.and  moment  is  very 
small.  The  segment  in  the  offset  position  has  the  largest  effect  on  the  force  and  moment 
coefficients  of  the  projectile.  Table  2  shows  these  aerodynamic  coefficients  for  the  segment  itself 
in  the  three  orientations.  The  axial  force  coefficient  is  increased  somewhat  when  the  segment  is 
offset  or  angled.  The  offset  segment  has  a  lai-ger  axial  force  than  that  obtained  for  the  angled 
segment.  As  for  the  normal  force  and  the  pitching  moment  coefficients,  the  effect  is  larger  with 
the  angled  segment  than  with  the  offset  segment. 
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Figure  10.  Surface  Pressure  on  Segments  (rear  view). 


Table  1 


Force  and  Moment  Data  for  a  Projectile 


Mach  =  2.5 

Alpha  =  0.0 

Segment  Position 

Centered 

Offset 

Ca 

+0.16 

+0.21 

+0.19 

Cn 

0.00 

-0.09 

+0.11 

CmD 

0.00 

+0.11 

-0.12 

Table  2 


Force  and  Moment  Data  for  a  Segment 


Mach  =  2.5 

Alpha  =  0.0 

Segment  Position 

Centered 

Offset 

Ca 

+0.23 

+0.22 

+0.25 

Cn 

0.00 

+0.04 

0.00 

CmD 

0.00 

-0.16 

0.00 

5.  CONCLUDING  REMARKS 

A  computational  study  has  been  imdertaken  to  compute  the  aerodynamics  on  a  multi-body 
projectile  system  consisting  of  an  ogive-cylinder  projectile  and  a  peg-shaped  trailing  segment. 
Flow  field  computations  have  been  performed  at  Mach  number  Moo  =  2.5  and  angle  of  attack  a  = 
0°  using  an  unsteady,  zonal  F3D  Navier-Stokes  code  and  the  chimera  composite  grid 
discretization  technique.  The  computed  results  show  the  qualitative  features  of  the  wake  flow 
field  for  the  projectile  alone,  as  well  as  with  the  segment  in  three  different  positions.  The 
predicted  flow  field  in  the  wake  regions  was  seen  to  change  dramatically  in  the  presence  of 
trailing  segment  in  the  wake  of  the  projectile,  especially  with  the  segment  in  the  offset  and  angled 
positions.  Computed  results  also  show  that  the  segment  in  the  offset  position  has  a  strong  effect 
on  the  flow  field  in  the  aft  region  of  the  projectile  upstream  from  the  base  comer  and  thus,  on  the 
aerodynamic  force  and  moment  coefficients  of  the  projectile.  The  aerodynamic  force  and  moment 
coefficients  of  the  segment  are  also  significantly  altered,  based  on  the  orientation  of  the  segment. 

The  computational  analyses  presented  in  this  report  show  the  capability  of  the  numerical 
technique  to  compute  the  aerodynamic  interference  flow  fields  associated  with  multi-body 
projectile  configurations.  The  numerical  technique  used  here  can  very  easily  be  extended  to 
determine  the  aerodynamics  of  multiple  segments  flying  in  the  wake  of  a  parent  projectile. 

Future  efforts  will  therefore  include  numerical  simulation  of  multi-body  projectile  configurations 
with  multiple  segments  at  (hypervelocity)  flight  speeds  of  interest. 
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